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Abstract 

We investigate the effect of static electron-phonon coupling, on real-time dynamics of spin and 
charge transport in vr-conjugated polyene chains. The polyene chain is modeled by the Pariser- 
Parr-Pople Hamiltonian with dimerized nearest-neighbor parameter to(l + ^) for short bonds and 
to(l — for long bonds, and long-range electron-electron interactions. We follow the time evolution 
of the spin and charge using time-dependent density matrix renormalization group technique, when 
a hole is injected at one end of the chain in its ground state. We find that spin and charge dynamics 
followed through spin and charge velocities, depend both on chain length and extent of dimerization, 
6. Analysis of the results requires focusing on physical quantities such as average spin and charge 
polarizations, particularly in the large dimerization limit. In the dimerization range 0.0 < 5 < 
0.15, spin-charge dynamics is found to have a well defined behavior, with spin-charge separation 
(measured as the ratio of charge velocity to spin velocity) as well as, the total amount of charge and 
spin transported in a given time, along the chain, decreasing as dimerization increases. However, 
in the range 0.3 < 5 < 0.5, it is observed that the dynamics of spin and charge transport becomes 
complicated. It is observed that for large 5 values, spin-charge separation is suppressed and the 
injected hole fails to travel the entire length of the chain. 



PACS numbers: 72.15.Nj, 72.80.Le, 71.10.Fd 



I. INTRODUCTION 



With vast advancements in technology, low- dimensional vr-conjugated organic systems 
in recent times, have found use in single-molecule electronic and spintronic devicesi""^ Un- 
til now, these materials have been used in devices such as organic light-emitting diodes 
(OLEDS) and organic thin-film transistors.-^-^ The vr-conjugated organic materials form an 
interesting class of strongly correlated systems in which there exists long-range electron corre- 
lations. Therefore, the low-energy physics of these systems is different from low-dimensional 
strongly correlated materials described by the Hubbard model. In order to propose and 
design vr-conjugated organic systems as components of electronic (spintronic) circuits, a 
proper theoretical understanding of the mechanism of charge and spin transport in these 
systems is thus necessary. Theoretical understanding of transport in many-particle systems 
with strong correlations requires appropriate techniques and formulations, mainly because 
transport is essentially an out- of- equilibrium phenomena. The advent of time- dependent 
density matrix renormalization group (td-DMRG) technique has vastly helped in addressing 
this issue.—"— Even so, spin and charge transport in vr-conjugated systems have not been 
addressed until recently using the td-DMRG method due to the fact that, most of the ex- 
isting td-DMRG algorithms are chiefly structured to handle short-range electron-electron 
interactions. Furthermore, those few that are capable of handling long-range interactions 
suffer from drawbacks of large computational resources for their study. The double time 
window targeting (DTWT) technique proposed by us^'', has been employed to address the 
issue of dynamics of spin and charge transport in vr-conjugated systems. 

The issue of coupling of the electronic and lattice degrees of freedom on the stability of 
the metallic state in one-dimensional conductors, has been of interest for a long time.— ~— 
One-dimensional systems such as polyenes, typified by trans-polyacetylene (t-PA), belong 
to the class of vr-conjugated molecular materials having linear (chain) topology. These 
materials are the simplest vr-conjugated systems that have been studied extensively both ex- 
perimentally and theoretically.— In these systems, electronic structure is strongly affected 
by electron-phonon interactions leading to dimerization, which is stabilized and enhanced 
by electron-electron correlations. As a result, these systems have dimerized ground state 
and their low-energy excitations are gapped, the extent of which depends both on the degree 



of dimerization and electron-electron interactions. Thus, it is to be expected that electron- 
phonon coupling will also influence the out- of- equilibrium dynamics of spin-charge separation 
in particular, and transport in general. Hence, an understanding of the role of dimerization 
on the dynamics of spin and charge transport in these systems is of considerable inter- 
est. Like most studies, we have also treated the electron-phonon coupling in the adiabatic 
limit (or Born-Oppenheimer approximation), i.e., the phonon degrees of freedom are consid- 
ered as slow, classical variables. In this limit, the sole effect of electron-lattice coupling in 
polyenes manifests as bond alternation, the extent of which is dictated by the alternation (or 
dimerization) parameter 6. A non-zero value of 6 implies that the nearest-neighbor hopping 
integrals alternate as to{l ± 6) . The object of this study is to investigate the effect of this 
bond alternation on the dynamics of spin and charge transport in one-dimensional systems 
typified by polyenes. The organization of the paper is as follows: In Sec. II, we discuss in 
detail, the model Hamiltonian and computational strategy used. Section III presents the 
results of our study along with discussions. In Sec. IV, we present our conclusions. 



II. MODEL AND COMPUTATIONAL METHODOLOGY 



The Pariser-Parr-Pople (PPP) model^"^ is appropriate for investigating the effect of 
dimerization on the dynamics of spin and charge transport in polyenes. In the second 
quantized representation, the PPP Hamiltonian reads as 

L-l 

^ppp = E - + 4i.<xc..) 

+ ^ Y^iifii - 1) + ^ Vij{hi - Zi){hj - Zj). 

i=l j>i 

Here, L denotes the number of carbon atoms in the polyene chain with open boundary 
condition, (Q,cr) creates (annihilates) an electron with spin orientation a on the i^^ 
carbon atom, to is the average transfer integral without dimerization and, < 5 < 1 is 
the bond alternation or dimerization parameter. The strength of on-site Coulomb repulsion 
between two electrons of opposite spins on site i is f/, and is the electron number operator 
for the same site. The term Vij represents the inter-site Coulomb repulsion between sites 
{i,j), with Zi being the on-site chemical potential of the i^^ carbon atom. Polyenes being 



homogeneous sp^ carbon systems, Ui = U aX all sites, and to maintain charge neutrality 
when a site is singly occupied, we also set Zi = 1 for all i. The inter-site interaction between 

2 

electrons on sites i and j, Vij, is interpolated between U for Vij = and f- for rij — oo 
by Ohno interpolation^ given by, 

Vij = 14.397(1.6348 + r^^.)-i/2. (2) 

where Vijs are in eV and r^jS are in A. Although there is an algebraic fall off in the inter- 
site potential Vij, all (equilibrium) properties such as optical gap, two photon gap and 
spin gap can be extrapolated to the thermodynamic limit. This is because the transfer 
part of the Hamiltonian is short-ranged and the effective correlation strength, Ves = {U 
- V12) is smaller than the band-width of the one-particle spectrum. Aside from this, the 
TT-coherence length (which is the length of the vr-system beyond which intensive properties 
such as excitation gaps, saturate) in such system is only about 20 sites long. Hence we 
expect that (non-equilibrium) properties such as spin and charge velocities, discussed in the 
subsequent section, also approach the thermodynamic limit. In this study we deal with 
polyene chains of 20, 30 and 40 carbon atoms, 6 is set to 0.0, 0.05, 0.07, 0.15, 0.3 and 
0.5, and the rest of the parameters assume standard values for the PPP model for t-PA 
and polyenes2i"^: = —2.4 eV, U = 11.26 eV, and 2n/3 bond angle between successive 
bonds. The PPP Hamiltonian possesses charge-conjugation and inversion symmetries, and 
also conserves total spin. Dimerization affects the transfer term and the distance- dependent 
electron-electron repulsions (Vij) only. It does not influence the on-site Coulomb repulsion 
(f/) between electrons. 

As 6 increases from 0.0 to 0.5, hopping integrals for partial double bonds get enhanced 
from —2.4 to —3.6 eV, while those for partial single bonds reduce from —2.4 to —1.2 
eV. As a result, the 1-norm of the Hamiltonian matrix increases, implying that = 
max[|ii^inax|, |-Emin|] also iucrcases. The dimensionless time-step a = Ec At of a numeri- 
cal scheme for solving the time-dependent Schrodinger equation, defines its stability region, 
with At being the time-step of evolutiort^. The value of a is constant for a given ordi- 
nary differential equation (ODE) solver, and hence as Ec increases. At has to be decreased. 
Thus, with increase in dimerization, the time-step for propagating the Schrodinger equation 
forward in time, has to be reduced for numerical stability. This is the scenario with com- 
monly used ODE solvers such as the Runge-Kutta (RK) schemes, the Crank-Nicholson (CN) 



method, and the multi-step differencing (MSD) techniques. Since the DTWT techniquei^ 
uses the MSD2 scheme^*^ for updating the Hilbert basis and the fourth-order RK technique 
for time evolution, with increase in dimerization the DTWT procedure becomes computa- 
tionally time consuming. Hence, we modified the DTWT algorithm by replacing these two 
time evolution methods with the Chebyshev-polynomial-based expansion of the time evolu- 
tion operator, U{At) = exp{—iH At) .—^ The Chebyshev-polynomial-based scheme has the 
advantage that the expansion of U (At) can be evaluated up to machine accuracy and is free 
from any time-step constraint. 

The Chebyshev-polynomial-based time evolution involves propagating the state | ipit)) by 
time-step At by approximating the discrete time evolution operator in terms of Chebyshev 
polynomials as 

p 

e-^^*^ J]a^T^(H), (3) 

m=0 

where, Tm(H) is the m**^ Chebyshev polynomial of the first kind, H represents the scaled 
Hamiltonian with eigenvalues ranging from [—1.0, 1.0], and the coefficients are given by 

a™ = (2 - 5^o)e-*^*^(-z)™J™(At/3), (4) 

where, 7 = (-Emax + -Emin)/2 and f3 = (Emax - -Emin)/2; -Emax, Emin are the maximum and 
minimum eigenvalues of the PPP Hamiltonian, and is the m*^ order Bessel function of 
the first kind. The necessity of scaling Hppp to H arises from the argument domain of the 
Chebyshev polynomials of the first kind T^ix); x G [—1,1]. The Chebyshev polynomials 
can be generated using the following recursion relation^^i^, 

m), (5) 

with the initial conditions, To(IH) \ ip) = \ ip), and Ti(]H) \ ip) = M \ ip) . However, since 
the coefficients am are known in advance [Eq. (4)], instead of using this forward recursion 
scheme, we use the "reverse" recursion algorithm proposed by Clenshaw^*^, which is more 
stable. The Clenshaw recursion requires P sparse-matrix vector multiplications (SMVMs) 
of the Hamiltonian H with the state vector | ip{t)). When P > 1 At (-Emax — -£'mm), the error 
decays almost exponentially'^-. In case of increase of 1-norm of the Hamiltonian matrix, one 
needs to either retain a higher P or reduce the magnitude of time step At. We increase the 
value of P with enhancement in dimerization. 



T^+i(H) I 



2HT^(H)-T„_i(H) 



The basic DTWT algorithm remains unchanged when ODE solvers such as the fourth- 
order RK and MSD2 schemes are replaced by the Chebyshev-polynomial-based expansion 
of the time-evolution operator. However, for a given number of retained density matrix 
eigenvectors (DMEVs) (m) and a given time step (At), results of the Chebyshev-polynomial- 
based DTWT td-DMRG algorithm are more accurate than the ODE-based version of the 
algorithm. This is because the ODE solvers have truncation errors associated with them; 
the truncation errors associated with the MSD2 and the fourth-order RK procedures are 
0[(At)^] and 0[(At)^], respectively.—'^'^ The Chebyshev-polynomial-based expansion of 
the time evolution operator, on the other hand, is free from such truncation errors and 
hence, can be evaluated up to machine accuracy by keeping P (number of SMVMs) > 
^At{Emax — Emin)-— However, the increase in accuracy is significant only at longer times; in 
the initial stages of time evolution, the accuracy is similar for the ODE and the Chebyshev- 
polynomial-based DTWT algorithms. 

The number of SMVMs associated with the basic MSD2 and fourth-order RK steps are 1 
and 4, respectively. For a system with sites, in the finite-system DMRG algorithm, for a 
single full sweep we have 4(A^/2 — 2) basic steps. Thus, total SMVMs for a single full sweep 
are 4{N/2 — 2) + 4. Since each single time window At is subdivided into p time slices of 
width At, for propagating the wave packet by a single window, the total number of SMVMs 
is 8p{N/2 — 2) + 4p. The first term corresponds to the MSD2 scheme for a full finite sweep 
of the site system over two time windows and the second corrresponds to propagation of 
the wave packet by a single time window. The total SMVMs per single time window are 
4p{N — 3)^. In case of the Chebyshev-polynomial method, the total number of SMVMs per 
single time window At is AP{N/2 — 2) or 2P(A^ — 4), P is dependent on At, and usually P > 
p. However, At in ODE methods should be very small (~ 10^^ fs) when the 1-norm of the 
Hamiltonian matrix is large, for reasonable accuracy, while in the Chebyshev-based method, 
At ~ 10~^ fs can be employed for similar accuracy. Thus, the number of time steps in ODE 
based methods (pode) is larger than in Chebyshev polynomial methods (pcp)- This leads 
to pcp X P ~ poDE and the two approaches require approximately the same computational 
time. 

In order to investigate the effect of dimerization on the dynamics of spin and charge 
transport, an up spin electron is annihilated from the first carbon site of the polyene chains 



of L sites in its half-filled ground state | thereby leading to an initial wave packet 

I m) 

|V^(O)) = ci,tl0y. (6) 

This wave packet is propagated in time by solving the time-dependent Schrodinger equation 
numerically, using the Chebyshev polynomial-based expansion of U{At). Using the time 
evolved wave packets | ipit)), site charge density [(nj(t))], and site spin density at 
time t are computed as, 

{fiiit)) = {^{t) I {ui,^ + I ^{t)), (7) 

(stit)) = I (n,,^ - n,,_.) I (8) 

These quantities give the dynamics of the injected hole in terms of its spin and charge 
degrees of freedom. The real-time dynamics of the initial wave packet is studied using a 
modified DTWT scheme, wherein adaptation of the Hilbert space as well as time evolution, 
are performed with the Chebyshev polynomial-based expansion of time evolution operator. 
The other parameters used in our study are: DMEVs retained, m = 300; time step for 
evolution At = 0.066 fs; total evolution time T = 33.0 fs. However, for the purpose of this 
study we focus only on the initial 15.0 fs, as this is adequate for our purpose. In the modified 
DTWT procedure, weights of the reduced density matrices of all the time-dependent wave 
packets are kept the same, unlike in the originally formulated scheme,— since we are using a 
large step size of 0.066 fs. Although time-dependent charge and spin densities are computed 
at all the sites, the quantities {ni(t)) and {hL{t)), and {sfit)) and are sufficient for 

investigating the effect of dimerization on the dynamics of the hole in terms of its spin and 
charge degrees of freedom; 1 and L correspond to the first and last sites of the vr-conjugated 
chain with L sites. 

III. RESULTS AND DISCUSSION 

The temporal variation in site charge density and site spin density for the terminal sites 
of a PPP chain of 40 sites are shown in Figs. 1 and 2, for various dimerization strengths 
6. The first significant dip in (fiLit)) and {sf^{t)) correspond to the times and r£, taken 
by the charge and spin of a hole injected at the first site, to reach the end of the chain, 
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FIG. 1: Temporal variation in charge densities, and (n4o(t)), at sites 1 and 40 in the PPP 

model, for different dimerizations, 6, for a chain of 40 sites. In each box, {ni{t)) is shown by thick 
black curve and (n4o(t)), by thin black curve. For 0.0 < 6 < 0.15, is indicated by arrow. 

respectively. The second significant dip in the time evolution plots of {ni{t)) and {sl{t)) 
occurs when the charge and spin degrees of freedom of the hole, get reflected back from the 
chain end and reach the injection site. These two dips correspond to times r^^. ^^^1 T2l, 
respectively. This is supported by the fact that r^^. are very close to twice the values 

of and r£, respectively. From these times, the charge velocity {'d'l) and spin velocity 



{i^D are calculated in a straight forward way: ^ 
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. From Figs. 1 and 2 it is 



observed that locating r^'^'^ is easy for dimerizations 6 up to 0.15. For dimerizations 6 = 0.3 
and 0.5, it is however very hard to locate these minima unambiguously. This may be due to 
the fact that the weak bonds, when very weak (large 6), do not easily transmit the charge 
or spin resulting in interference of the wave traveling forward with the reflected wave. This 
effect is seen more in the longer chains as the charge or spin needs to travel through many 
weak bonds. 

Normally, charge and spin velocities are discussed in the bulk limit. From our earlier 
studies on the one-dimensional Hubbard model, we have noticed that it is not possible 
to relate the analytically obtained charge and spin velocities "^p/a, with 'i9^/2L obtained for 
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FIG. 2: Temporal variation in spin densities, (sf(i)) and (s|o(*))) at sites 1 and 40, in the PPP 
model, for different dimerizations, 5, for a 40-site chain. In each box, (sf(t)) is shown by thick 
black curve and (s4Q(i)), by thin black curve. For 0.0 < 6 < 0.15, t£ is indicated by arrow. 

finite small systems. This is because, the analytic velocities are obtained from the exact 
excitation spectrum as group velocities in the limit of momentum g*— )■ 0,^— while in our case 
pertain to a state which is not an eigenstate of the Hamiltonian governing the dynamics 
of the system. Furthermore, dimerization introduces a gap in the excitation spectrum and 
there is no low energy theory which completely separates the spin and charge degrees of 
freedom. This is also true for the one- dimensional PPP model [Eq. 1] for which however, 
no analytical expressions for the spin and charge velocities are available in the literature. 
Hence, unlike i!^p/a, ^^^^ depends on system size. This dependence is more pronounced in 
the Hiickel and Hubbard models, compared to the PPP model, due to the absence of long- 
range interactions.— The ratio of the velocities, due to finite size effects is system size 
dependent. We need to go to longer chains (longer L) for this ratio to be size independent. 
For the PPP model however, this dependence is weak (see Fig. 3). 

In Table I we present the times taken by the charge and spin of the hole to propagate from 
one end of the chain to the other end, as well as, the charge and spin velocities, and their 
ratios. The times taken by the spin and charge to reach the end of the chain approximately 




scale with chain length. However, the velocities of both charge and spin increase slightly with 
increasing chain length, which is due to weak finite-size effects. In Table I, we also present 
the spin and charge velocities extrapolated to infinite system size for different dimerization 
strengths (5 = 0.0, 0.05, 0.07 and 0.15), as well as, the ratio ^^/^l^- We see that the 
ratio '&'^/'&lo decreases with an increase in 6. In Fig. 3 is shown the dependence of the 
ratio of charge to spin velocities as a function of dimerization, for the chain lengths studied. 
We find from the plot that the ratio decreases as 6 increases. Indeed, within the error of 

h is 

resolution of , even for 5 = 0.3 and 0.5, it appears that this trend continues and the 
ratio of the velocities approaches the non-interacting value of 1.0. This can be contrasted 
with the push-pull polyene systems studied by us,— where the ratio of the spin and charge 
velocities is independent of the strength of the push-pull groups. The above feature, namely, 
dependence of 'd'l/'d]^ ^) cannot be attributed to change in interactions brought about 
by geometry changes as a consequence of increased dimerization. This is because these 
changes are small and it has been shown in earlier studies^'^ that the contributions to the 
energy gaps between states due to changes in interaction parameters, V^j, accompanying 
small bond length changes (due to change in 5), is rather small. Therefore it appears that 
in dimerized models, the change in transfer integrals has a stronger role to play than the 
long-range interactions. 



In order to investigate in detail the issue of the ratio ^ tending to the non-interacting 
value of 1.0 as 5 — )■ 1.0, we have focused on the time dependence of the charge (spin) po- 
larization defined as, Pc(s)(^)- These are calculated from the normalized site charge (spin) 
density pf'^\t), with background correction. Employing these renormalized quantities sup- 
presses the large fiuctuations seen in {riiit)) and {slit)) allowing one to focus on the essential 
behavior of spin and charge transport. The quantities plit) and plit) are defined as 

... ^!^,M^c_. {sm)-s 

where C and S are the average background charge and spin densities respectively, and are 
given by, 

^ = ^ 1 % I 0°.); s = ^ J2(<p% I s^^ I <!>%). (10) 
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FIG. 3: Variation in the ratio of charge and spin velocities {'&^/'&%) with dimerization 0.0 < 5 < 
0.15, in the PPP model, for different chain lengths. 



For a half-filled ground state belonging to the S^^^ = subspace, C = 1.0 and S = 0.0. Using 
these observables, we define time-dependent charge and spin polarizations as, Pc(s)it) = 
p^j'^\t)rj. The time-dependence of the charge and spin polarizations qualitatively reflect 
the center-of-mass movement of the charge and spin peaks.— However, these polarizations 
are gauge dependent as the system is having a net charge. Hence we have computed these 
by choosing as origin, the central site of a site chain (A^ = 20, 30, 40). 

In Fig. 4, we have plotted the variations in charge and spin polarizations with time, for 
a polyene chain of 40 sites. Motion of the charge (spin) across the two ends of the chain 
leads to oscillations in the temporal variation of Pc[s){i)- The larger the amplitude of the 
oscillations, the higher is the probability of charge (spin) traveling between the chain ends. 
It is observed from both the curves that, with an increase in dimerization, the oscillation 
amplitude decreases. This indicates that electron-phonon coupling reduces the probability 
of both the charge and spin degrees of freedom to travel between the two ends of the chain. 
This observation indeed supports our earlier observation that it is difficult to locate tI^^ for 5 
= 0.3 and 0.5 merely because, weak bonds fail to transmit the charge or spin easily. In other 
words, the charge and spin for large values of 6, might not even reach the opposite end of the 



TABLE I: Variation in the times r]j, t£ (fs), and velocities ■i?^, ■i?!^ {A/ is), and ratio of velocity of 
charge to velocity of spin ■d^/'&j^, in the PPP model with dimerization, for different chain lengths, 
L = {N — 1)7-0 - (5; AT is the number of sites and ro = 1.397 Ais the C = C bond length in a regular 
polyene chain. 
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chain, especially for long chains. We also notice that an increase in dimerization suppresses 
spin-charge separation, which is manifested by the fact that Pdt) and Psit) closely follow 
each other. Thus, indeed the ratio of-^^— >-1.0as(5^1.0. In other words, with an increase 
in electron-phonon coupling, spin-charge separation is suppressed in the PPP model. 

However, the question about fate of the charge and spin of the injected hole for large 




5 10 
time (fs) 



5 10 
time (fs) 



FIG. 4: Temporal variation in Pc{t) (dark curves) and Ps{t) (light curves) for a 40 site PPP chain, 
for different dimerizations, 5. 



values of 6, still remains. In order to investigate this issue, we have computed 

dNnLjt) _ N^{t)-N^{t) 
dt At 



dt 



At 



(11) 
(12) 



Here, Nn/L{t) = EjeR/LinM and 5^/Ji) = E,e«/L(«|(^)); R e [N/2 + 1,N] and L e 
[1, N/2\. These quantities provide information about the amount of charge and spin of the 
hole that is transported from the left-half to the right-half of the chain. In the initial state of 
the system, both the charge and spin of the hole reside solely in the left-half of the system. 



Hence, 



and 



at time t = have the values 1.0 and 0.5, respectively. With time. 



dt dt 

as the charge and spin travel from the first to the last site of the chain, these quantities 
change sign and depending on the strength of electron correlations, attain values close to 
— 1.0 and —0.5. When the charge and spin travel back from the last to the first site, these 
quantities again change sign and reach values close to -1-1.0 and -1-0.5. This thus indicates 
that, '^^^^^^^ and oscillate between ±1.0 and ±0.5, respectively. The amphtude of 

oscillation of these quantities with time gives a measure of the amount of charge and spin 
that has been transported from the left-half to the right-half of the system. If under some 
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FIG. 5: Variation in '^^^^^^ (dark curves) and '^^^^^^ (light curves) with time, for a 40 site PPP 
chain, for different dimerization strengths: 6 = 0.0, 0.05, 0.07, 0.15, 0.3, and 0.5. 



circumstance the charge (spin) fails to travel between the ends of the chain, time evolution 
of these two quantities does not not show any sign change. 

Figure 5 depicts the time evolution of '^^^^^^^ and ^^^^^^ for a polyene chain of 40 sites, 
for different dimerization strengths. It is observed that these quantities indeed oscillate 
between ±1.0 and ±0.5, although the amplitude of oscillations decreases with time, as 5 
increases. For small dimerizations {S = 0.05 and 0.07), it is observed that the temporal 
variations in '^'^^^^^^ and '^^^ are similar to those in the undimerized case, and the number 
of times "^^^^^^^ and '''^^^^^ changes sign is the same for all the three cases. However, the 
number of "sign changes" in the time propagation of "^^^^^^ is more than that of — 
indicating that the charge degree of freedom moves faster compared to the spin degree of 
freedom which is merely a manifestation of spin-charge separation. Thus, for small values of 
clcctron-phonon coupling the motion of the charge and spin of the hole do not get hindered 
due to interferences. 

However for 5 — 0.15, 0.3 and 0.5, the time evolution profiles of "^^^^^^^ and '^^^^^^ are 
completely different from those obtained with small values of S. It is observed that the 
amplitude of oscillations diminishes with time indicating that the amount of charge and 



spin transported across the chain decreases as electron-phonon couphng increases. This 
supports our observation that an increase in 6 leads to a decrease in the height of the 
minima, in the time evolution profiles of {nL{t)) and (s|^(t)). Also, the number of times 
that ^^^^ and chan ges sign decreases, which reflects that strong electron-phonon 

coupling significantly reduces the probability of to-and-fro motion of the charge and spin 

dS^ (t) 

along the chain. For 6 = 0.3 and 0.5, however, — does not change sign within the time 
scale of our studies. This depicts that, for these values of 6, only the charge degree of the 
hole can travel across the length of the chain, albeit slowly and in a very small amount. The 
spin of the hole fails to reach the end of the chain and "gets trapped" within the left-half of 
the chain. This is the reason why we find that locating r^*^* for large values of 6 is difficult. 

Bosonization studies on one- dimensional spinful Hubbard and Luttinger liquid models, 
at half-filling^""—, indicate that 6^0 introduces a backward-scattering operator into the 
bosonized Hamiltonain, which is associated with momentum transfer q ~ 4k p. Because of 
this umklapp scattering term, as dimerization increases, propagation of the charge and spin 
degrees of freedom of the injected hole between the ends of the chain become hindered, and 
for large values of 6 (0.3 and 0.5), fail to reach the chain end (L). Furthermore, periodic 
perturbation of the lattice such as Peierls distortion, which introduces a gap in the excitation 
spectrum, has been shown to couple the charge and spin sectors of the one-dimensional 
(bosonized) Hubbard Hamiltonian.-^ These two results, obtained with respect to spinful 
fermion models with short-range electron-electron interactions, seem to support both our 
observations in the context of the PPP model, namely, an increase in dimerization suppresses 
spin-charge decoupling, and inhibits the oscillatory to-and-fro motion of the charge and spin 
of the injected hole. 



IV. CONCLUSION 



To conclude, we have shown that dimerization modifies the dynamics of spin and charge 
transport in the PPP model. For small dimerizations the system shows a smooth decrease 
in the ratio of charge to spin velocities, although both the velocities decrease slightly. For 
large values of 6, the situation is different. It is difficult to obtain the charge and spin 
velocities as strong interference effects smear out their motion. However, a more careful 



analysis based on net charge transport across the chain shows that both the charge and spin 
stay almost localized for the duration of our time evolution. Thus, although dimerization of 
polyenes predominantly changes the one-electron part of the Hamiltonian, it strongly affects 
spin-charge separation. Large values of dimerization reduce the amount of charge and spin 
transported across the system as well as, suppress spin-charge separation, as evidenced by 
the decrease in ratio and by closely following time evolution profiles of spin and charge 

polarizations. 
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